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1. Introduction 

Elementary reactive process between two molecules A and B usually begins when their 
relative positions satisfy certain relations. In the simplest case of two molecules with 
spherical symmetry, one may consider that the reactive process starts when the relative 
distance of their centres is equal to the range of their reactive potential. This is the 
case for instance in the celebrated model introduced by Smoluchowski P and improved 
by many authors |2j. In a more general model, for given values of the coordinates 
characterizing A, the reaction can only take place if the coordinates of B are located 
inside a reactive region T(A) which, obviously, can be fairly complicated. 

The probability that B has not reacted up to time t is clearly connected to the 
residence time of B inside T(A) up to time t. This fact seems to have been first recognized 
by Wilemski an Fixman jH]. Indeed, in their seminal work, the inverse reaction rate 
constant for the A + B — > B diffusion limited reactions is given by the mean residence 
time (RT) of the particle B inside T(A). We address the reader for more details to 
the works by Doi [3], Temkin and Yakobson [3], Zhou [S| and Rice 0. The residence 
time inside a finite domain S is defined as the total time T(t) spent by the Brownian 
particle up to time t within S: contrary to the first exit time of S, the RT accounts for 
an arbitrary number of exits from and enters into S. The properties of the RT for an 
individual Brownian particle have been extensively studied [3 El El EB CEH EE21 CE31 IE] • 

The purpose of this paper is to develop the notion of residence time in the case 
of a reaction between two particles, activated by immobile catalytic sites. For such 
reactions, A + B + C — > B + C, where the specie C stands for the catalytic sites, the 
calculation of the effective reaction rate in a Wilemski-Fixman kind approach actually 
involves several functionals of time which can be identified as RT of Brownian particles 
[To] . One of them is the mean time spent simultaneously by the two reactants A and 
B inside a reactive region centered on a catalytic site C, for which the RT is averaged 
over all initial positions of A and B inside the reactive region. 

Recently (Ej, this kind of joint residence time properties has been studied for an 
arbitrary number N of Brownian particles, but only in the two-dimensional situation 
and for very particular initial conditions. Here, we focus on the case of two particles 
(N = 2). First, it will be shown that the mean residence time can be calculated exactly 
for general initial conditions in the two dimensional case. Then, in order to obtain an 
explicit reaction rate for the reaction A + B + C-+B + C in 3-d, we will also study 
the three dimensional situation when particles start inside the reactive region. We will 
show that the mean joint residence time of A and B on a site C can be calculated 
explicitly in several situation of physical interest. It is remarkable, in particular, that 
this problem can be solved exactly in 3 dimensions for arbitrary initial positions of the 
reacting particles. 
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2. The model 

We consider 2 independent Brownian particles diffusing mad dimensional space, with 
d > 2, and we assume that they can only react on an immobile catalytic site, which 
will be represented by the sphere D(0,R) of center O and of radius R. We note Tj(t) 
the positions of these particles at time t (j = 1,2), the initial positions, and Dj 
the corresponding diffusion coefficients. We introduce I5 the indicator function of the 
sphere D(0,R), that is 



Mr) 



1 if re D(0,R) 
otherwise 



(1) 



(2) 



The joint occupation time of the sphere D(0, R) is defined by 

T( t ) = f&t> fn 15MO) 

Thus, the mean joint occupation time at time t is 

tit) = ((T(t))) = f (n<lsMO)>;) < (3) 



that is 



m = 1 n [j *.i*jp^^ j j dt ' (4) 

From now on, we will limit ourselves to the long time limiting value of /i(t), which 
will be simply denoted as /x 

1 r°° ^ If. I II™. — r ■ I 



» " fc ^ = (i6^D^ r d ^ 2 n / dr i^( r i) 



(IGttW.D^ Jo ™ n^^ron;-^ , ^ 

where we did the change of variables y — 1/t' in the integral. 

In the Wilemski-Fixman approach of the bimolecular reactions A + B — > B [HHUE], 
the quantity which is actually involved in the expression of the kinetic constant is the 
mean RT (/j) of a Active particle inside a reactive sphere of radius + vb, averaged 
over the initial position of the particle inside the region, diffusing with D = D a + -Db- 
More precisely, the overall rate constant k is given by the inverse addition law 
111 
k ko kA 

where hp is the contribution due to diffusion of particles and k^ is associated to the 
reactive process itself. It is shown that 

kD = w (7) 

where V is the reaction zone volume. 
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Note that this Wilemski-Fixman approach actually appears as the leading order 
of a perturbative theory with respect to which, as an unperturbed motion, takes 
the propagation of independent particles^Zj. In the framework of this approximation, 
the mean residence time (//) is calculated by considering that the particles perform 
independent Brownian motions. 

Our main purpose here is to compute the mean value (//) for the present catalytically 
activated reaction A + B + C^B + C, which now involves two diffusing particles and 
a disjoint reactive region. Therefore, (/i) will now be called the mean joint residence 
time inside a region. The sequel of the paper is structured as follows. In section El 
we present results concerning the case when all the particles are initially at the origin 
O. These expressions are valid for an arbitrary dimension of space d. Next, we study 
more general initial conditions in dimensions d = 2 and d = 3 which are the physically 
relevant dimensions of the problem. In section El we discuss the d = 2 situation in two 
special cases: firstly, when all the particles are at t = on the boundary of the sphere 
D(0, R). Secondly, when one of the two particles is initially at the origin while the other 
one can be everywhere. In section El we solve the d = 3 situation in full generality, that 
is for arbitrary initial conditions of the diffusing particles. Therefore, this result gives 
the mean joint residence time involved in the bimolecular reaction activated by catalytic 
sites, whereas the other initial conditions mentioned previously can be related to specific 
reactive processes. 



3. d-dimensional case when the two particles are initially at the origin 

Physically, this situation corresponds, for example, to a dissociation recombination 
catalytically induced reaction. 

We start by integrating Eq. © over y, which leads, after rescaling of the spatial 
integration variables to 

where 

I(D 1 ,D 2 ) = dx dy 2 ," 1 (9) 

Jo Jo (x z + y l ) a 1 

and ad is the area of the d-dimensional hypersphere of unit radius. 

We now proceed to the calculation of I(Di,D 2 ). First, we write I(Di,D 2 ) = 
J(D U D 2 ) + J(D 2 ,D 1 ) with 

J(D U D 2 ) = / dxdy % (10) 
Jv (x z + y l ) a 1 

and V is the domain V = {x G [0, ^=], y G [0, ^=], y < \f%x}- 

Changing to the new variables z — x, t = y/x in J(Di, D 2 ), we obtain 



J(D h D 2 ) = I dzz dt- 



t 2 ) 



d-l 
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1 ( Di\ d / 2 /, d 1 d D x \ 

2 Fx d- !,-;! + -; --i , (11) 



8d£>i V-D 2 / V 2 2' D> 2 

2^1 being the standard Gauss hypergeometric function. 
Thus, 



a 2 d R 2 T(d- 1) f 1 ADA 1 (D 



A*= * o,_, 7 7T^ 7T +-R-g hr (12) 



8d7r d iDi \£) 2 / D 2 \D 1 
where 

<7(:r) = x rf / 2 2 Fi (d - 1, ^; 1 + ^; (13) 
Using the expression cr^ = d7r rf / 2 /T(l + d/2) as well as the duplication formula 

r(i _ 1) = ^rf^iV^ (") 



> V 2 / \2 

we finally find 



,4 2 «iffi(> §H (1B) 



7T r(f + i) I A UV A VA 

Specializing to the dimensions d = 2 and d = 3, we get the following explicit 
expression of the mean residence time fi as a function of i?, D\ and A : 

• for d = 2 : g(x) = ln(l + x), thus 

^?l^lnfl + ^ + 1 Llnfl + #n (16) 



4 [D 1 V IV D> 2 V Di 

This result has already been obtained in ref fH] by a slightly different approach 
3 

for d = 3 : g(x) = - arctanf-v/x) — , v - , and 




' 1 + x 

R 2 ( 1 / [d[\ 1 
u = — < — arctan W — H arctan \ — . > 17 

We now briefly comment on the expressions (fTSJl and (fTTj) in the special case when 
the particle 1 moves much faster than the particle 2, that is when D\ S> D 2 . 
In that case, the two dimensional expression Eq. (jl6j) becomes 

R 2 , fDi\ 

^~^77rM "FT J ■ (18) 



The form of this last result appears, from the physical point of view, to be very 

reasonable. As a first approximation, we can say that this joint occupation time is 

given by the joint occupation time before the first exit time of the slow particle out of 

R 2 

the disc. This time is of order t 2 = — — . Meanwhile, the fast particle leaves the disc and 

D 2 

comes back several times. The time spent at each return of the particle 1 inside the disc 
R 2 

is of order t\ = — — , and the order of magnitude of the number of returns at time t 2 of 
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't_2 
A. 



the particle one inside the disc is given by In ( — ) , in the limit of long t 2 times. As a 



consequence, we expect 

"~^» ( ' in (!)~iHt) < i9) 

which reproduces the dependence of \x as a function of the diffusion coefficients D\ and 
D 2 given by Eq. |T6|) . 

Concerning the three dimensional expression Eq. (jl7)l . the limit D\ ^> D 2 leads to 

" wi- < 20 > 

This is clearly the expected result: knowing that the three dimensional Brownian motion 
is transient and that the particle 2 leaves the sphere after a very long time, the mean 
joint residence time identifies itself to the residence time of the only particle 1, which is 
precisely given by 

We now examine more general initial conditions, and we restrict ourselves to the 
physically important cases corresponding to the two-dimensional and three-dimensional 
situations. 



4. Two-dimensional situation 



Integrating first with respect to the angular variables in Eq. (0) 

" = Ilk f ^ jS (/„ V 3 'o (0) exp (-^ + (rf )»))) 

where Io stands for a modified Bessel function of zeroth order. We have not been able 
to compute explicitly this integral for arbitrary initial conditions, but only in special 
cases: the first one is the case when the two particles are initially on the boundary of 
the disc and the second one the case when one of the two particles is initially at the 
origin. 



4-1. Case when the two particles are initially on the boundary of the disc 
In that special case, the spatial integrations can be done explicitly. 

os ( vR 2 \ /', t (yR 2 v.j\ ( w 2 , \ 

- K exp ( —77^- ) / dUjUjlo [ n ^ ] exp 1 



ADj) Jo J J V 2D j J V AD 



j . 



where we have used the fact that 



f 1 dxxe~ ax2 l (2ax) = f e a - e~ a I (2a)) (23) 
Jo 4a v ' 
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Finally, 



R 2 y 



'yR 2 



4 Jo y 2 



(24) 



7=1 \ \ "^3 J j / 

In the special case when the two particles have the same diffusion coefficient 
D\ = D2 = D, this equation becomes 

R 2 r° dt / . . x 2 



a* r°° at / tr , , 
p 8L> Jo t 2 V ow 



where the last integral can be numerically evaluated as 

r°° dt / , . \ 2 

y ^ 



1 - e-*I 



1.06 



(25) 



(26) 



Case w/ien one of the two particles is initially at the origin 

We consider now a more general situation corresponding to the case when one the 
two particles is initially at the origin while the other one can be everywhere. This 
situation can be a simple model for a third order reaction between an immobile site C 
and particles 1 and 2, if C and particle 1 initially constitute a complex molecule which 
should be dissociated to react with a chemical species 2 (provided that the distances 
from 1 and 2 to C are less than R). 

The residence time of 1 and 2 within the catalytic region is then 

1 



^ = ^r^f dririexp 



dr 2 r 2 I 



(0) 

yr\ V 2 



cxp 



2D 2 ] "V 4L> 2 
We perform first the integral with respect to y: 




(27) 



dy exp 



-y 



AD, 



In 



yrf ] r 2 
2Do 



(oh 2 N 



4Di 



(r 2 r 



(oh 2 
2 J 



AD\ 



Performing next the integration over r 2 and lastly over r±, introducing the dimensionless 

r (o) 

variables x = and w = j^, we find after straightforward but lengthy calculations 
that 



if x < 1 



if x > 1 



/* = ^ {^-(1 - In A 2 ) + i-(l - x 2 - In AO - ±- ((1 + * 2 + «;) 2 - 4* 2 ) 1/2 } (29) 



^ = T {zT (1 " HA2) + zT (a;2 " 1 " ln(xMl)) " >f ( (1 + x " + w)2 " 4x2 ) V2 } 
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Ay = |(1 + x 2 + wf - 2x 2 - (1 + x 2 + w) v /(1 + x 2 + m;) 2 -4x 2 | (31) 



and 

A 2 = — |(1 - x 2 + w) 2 + 2wx 2 - (1 - x 2 + w) v /(l + x 2 + w) 2 -4x 2 | (32) 

In the special case when the two particles have the same diffusion coefficient 
Di = D 2 = D, this last expressions simplify in 

A* = ^/(*)> (33) 



where 



2 / — i — I — 7 i /a; 4 + 8 + 4 v / x T T4\ 
x l - V x 4 + 4 + In I I 1 < x 



2 - x 2 - Vx 4 + 4 + In (x 4 + 8 + 4Vx 4 + 4) x < 1 

(34) 

In the spirit of formula (|IJ), it is interesting to obtain the mean value of /i with 
respect to the initial position of r 2 . Assuming that the particle 2 is initially uniformly 
distributed in the disc D(0,R), we get 



tt/? 2 
/IV 

16 l/V" \DJ ' V-Di 



+ (35) 



where 



&i(x) = a(x) — x + 2 In 
6 2 (x) = -2a(x) + 2 In 



2(x + a(x)) 



x(x + 2 — a(x))(— x + a(x)) 
2 \ 



^x 2 + 4x + 2 - (x + 2)a(x) / 
a(x) = yj(x{4 + x)) (36) 
In the case of equal diffusion coefficients, it is given by the very simple formula 

(?) « 0.28 ^. (37) 

-R 2 / D \ 

In the limit D\ 3> _D 2 , we obtain (/i) ~ In ( — I, this result admitting the 

same interpretation than the one given after Eq. (jl6j) . 



5. Three dimensional situation 



We now turn to the most challenging situation, that is the three dimensional case for 
arbitrary initial conditions. 
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Integrating first over y in Eq. © and denoting for simplicity = v and r 



2 



U. 



we have 



1 2 



n dr ' l5 < r <>— 

V 4D, / 



(27r) 2 /•*, /•« . /•«•,. r 2 ?^ sin 6>i sin 6> 2 

dri / dr 2 / dfc>i / d# 2 - 



(167rAA) 3/2 Jo ./o ./o ./o ^( r 2 + v 2_ 2ri? ;cos6l 1 ) + ^(rl+u 2 -2r 2 ucose 2 )^ 

(3J 

It turns out that all the integrations can be performed exactly, and /i can be 
expressed in terms of elementary functions of the initial distances v and u of particles 
1 and 2 to the center of the catalytic site. These (rather lengthy) expressions are listed 
in the Appendix. 

Averaging over the initial positions of particles (which are supposed to be uniformly 
distributed in the reactive zone), we find 

(4tt) 2 J*dvv 2 Jo* dm/VM 



(4tt#73) 2 

R2 { ^m + ^m(^)\ (39) 



20vr IDi \D 2 J D 2 \D 1 



where 



2-101n(l + x) ln(l + x) , . — , 7 1/9 
77^ - 2 1 J 9 ; + 16 arctan(v^) - -x 1/2 (40) 



2-1/2 x 3/2 



Note that in the special case when one of the two particles, say 2, is immobile 
(-D2 = 0), following the Wilemski-Fixman jHl lU Ej approach and Eq.fjZJ), we obtain the 
expression 

k D = |(47rZ)ii2) (41) 
o 

given in the paper of Do'i 0, which differs from the Smoluchowski constant only by the 
factor 5/6. 

Note also that in the special case when the particle 2 is immobile and initially within 
the catalytic region, we recover the expressions given in ^2] fo r the mean residence time 
of a single Brownian particle within a sphere. 



6. Conclusion 



We have computed the mean joint residence time of two particles in a sphere, for very 
general initial conditions. In particular, we have obtained in all cases the dependence 
of this residence time as a function of the diffusion coefficients of the two particles. 
Such results can be especially useful for describing kinetics of a bimolecular reaction 
activated by catalytic sites. For this reaction in three dimensions, the effective reaction 
rate actually involves the mean joint residence time of two reactive particles inside a 
region centered on a catalytic site. This mean time, being averaged over all initial 
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positions of the two reactive particles inside the reactive region, is explicitly given by 
Eq.(|3I2). This is our main result, which appears as a generalization of the expression of 
the reaction constant given by Wilemski and Fixman for a bimolecular reaction. 

7. Appendix : Expressions of \i in three dimensions 



Considering the expression ()38l of /i, we integrate over angular variables and obtain : 

r\r\ sin 9\ sin 9 2 D\D 2 r\r 2 



d0i / d9 2 



Di 



r\ + v 2 — 2r 1 vcos9i) + j^(r 2 + u 2 — 2r 2 ucos9 



Avu 



In 



{D 2 {r 2 + v 2 ) + Di(r\ + u 2 ) + 2D 2 vr l - 2D 1 ur 2 )(D 2 (r 2 + v 2 ) + D x (r\ + u 2 ) - 2D 2 vr x + 2D lU r 2 ) 
(D 2 (r{ + u 2 ) + Di(r% + u 2 ) - 2D 2 vr l - 2D 1 ur 2 )(D 2 (r 2 + v 2 ) + D x {r 2 2 + u 2 ) + 2D 2 vr l + 2D x ur 2 ) 

(42) 

Integrating next by parts over the radial variables r\ and r 2 , we finally find for initial 
positions v and u 



where 



fj.(v,u) 



1 



96iruv 



/ 1 



+ 8m; ( — ft(u, w; £>i, L> 2 ) + ttM". «; D 2 , Di) + 24iwi? 2 — «; A, D a ) + 7^0, u; D a , A) ) + 
L)\ L> 2 J \-L> 2 D\ 



Dju 4 + Dy - 6D 1 D 2 n 2 v 2 - QR 2 (D\u 2 + D 2 v 2 - D 1 D 2 (u 2 + v 2 )) - 3R l (D 1 + D 2 f V ' Dl,D ' 2 ' 



(AA) 3 / 2 



+ 



+ 



8R 3 



{D.D.f/ 2 



D 2 u9(u, v; D x , D 2 ) + D 2 2 v9(v, u; D 2 , + v {u, v; D u D 2 



(43) 



f(u,v;Di,D 2 ) = arctan 



arctan 



h(u, v; Di, D 2 ) = arctan 



+ arctan 



ip(u,v;D 1 ,D 2 ) 



lh(R 



£>i (R 
^(R 



D 1 {R-u 
H^(R + u 



D 2 (R 
^(R 



D 2 (R 



arctan 



arctan 



arctan 



arctan 



D! (R 
^(R 



u 



D 1 (R + u) 

flh (R + u) 
D 2 (R + v) 

FdI(r-u) 

D~ 2 (R + v) 



(44) 



+ 



(45) 



In 



((L>i + D 2 )R 2 + D lU 2 + D 2 v 2 - 2D lU R + 2D 2 vR)((D 1 + D 2 )R 2 + D lU 2 + D 2 v 2 + 2D lU R - 2D 2 vR)\ 
((A + D 2 )R 2 + D x u 2 + D 2 v 2 - 2D 1 uR - 2D 2 vR)((D 1 + D 2 )R 2 + D lU 2 + D 2 v 2 + 2D lU R + 2D 2 vR) J 

(46) 



9(u,v;D h D 2 ) = 

((A + D 2 )R 2 + D lU 2 + D 2 v 2 - 2D lU R + 2D 2 vR){(D 1 + D 2 )R 2 + D x u 2 + D 2 v 2 + 2D lU R + 2D 2 vR) 
((A + D 2 )R 2 + D x u 2 + D 2 v 2 + 2D lU R - 2D 2 vR)({D l + D 2 )R 2 + D lU 2 + D 2 v 2 - 2D lU R - 2D 2 vR) 



In 



(47) 
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R 2 (— + — ^ u < R and v < R 

R 3 „ 

R < u and v < R 



r](u,v;D 1 ,D 2 ) 



1/3 
1/3 



uD 2 
R 3 



u < R and R < v 
otherwise 
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